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An upper limit on the total annihilation cross section of dark matter (DM) has recently been 
derived from the observed atmospheric neutrino background. We show that comparable bounds 
are obtained for DM masses around the TeV scale by observations of the diffuse gamma-ray flux 
by EGRET, because electroweak bremsstrahlung leads to non-negligible electromagnetic branching 
ratios, even if DM particles only couple to neutrinos at tree level. A better mapping and the partial 
resolution of the diffuse gamma-ray background into astrophysical sources by the GLAST satellite 
will improve this bound in the near future. 
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I. INTRODUCTION 

One promising way to detect dark matter (DM) 
is indirectly via its annihilation (or decay) prod- 
ucts. The DM annihilation productsbarring ex- 
otic models with additional stable and relatively 
light particlesare standard model (SM) particles, 
although with model-dependent branching ratios. 
Using atmospheric neutrino data, the authors of 
Ref. [l| derived an observational upper bound 
on the annihilation cross section ((Tannic) of any 
DM candidate, assuming that it annihilates only 
into the least detectable final states in the SM, 
namely neutrinos. Allowing only couplings to neu- 
trinos might be not only a conservative assumption 
needed to derive this bound, but could be real- 
ized in nature: Possible DM candidates, like the 
Majoron, with this property exist. Moreover, the 
bound on the diffuse gamma-ray background from 
EGRET observations [3, H, |3| translates into ex- 
tremely restrictive limits on the branching ratios 
in electromagnetic and hadronic DM annihilation 
channels. Therefore, models with high annihila- 
tion rates proposed to solve the cusp problem of 
conventional cold DM (see e.g. [H]) and DM masses 
rax above 0{GeV) are likely to require either fine- 
tuning or should couple the DM particle only to 
neutrinos. 

The latter possibility has already been invoked 
in exotic scenarios explaining the origin of ultra- 
high energy cosmic rays. Reference [^proposed 
that supermassive relic particles decay only into 
neutrinos, thereby contributing to the ultrahigh 
energy cosmic ray flux through the Z burst mech- 
anism and avoiding at the same time constraints 
from the diffuse gamma-ray background. However, 
the authors of Ref. 0] showed that electroweak jet 
cascading leads to a non-negligible electromagnetic 
branching ratio and rules out these models. 

In this work, we extend this argument to anni- 
hilating dark matter of lower mass, showing that 
this mechanism combined with the limit on the 



diffuse gamma radiation by the EGRET satellite 
provides competitive observational constraints on 
Wannv) for masscs around the TeV scale. Future 
observations of the diffuse gamma-ray flux by the 
GLAST satellite should improve these bounds. We 
also find that the strongest and most robust way 
to constrain ((lannw) is to use the DM signal as- 
sociated with the galactic halo, instead of the dif- 
fuse flux from cosmologically distributed dark mat- 
ter. We comment on the possibility to improve the 
neutrino bounds as well by exploiting the strongly 
peaked angular distribution expected from annihi- 
lations in the galactic dark-matter halo. In Sec. HIl 
we discuss the properties of dark matter relevant 
here, while Sec. IIIII is devoted to the data used 
to derive the bound. The bound is presented and 
commented upon in Sec. [iVl In Sec. |Vl we discuss 
possible improvements and flnally conclude. 



II. THE DARK MATTER INPUT 

In Ref. fH, the expected dominating contribu- 
tion to the diffuse neutrino flux was estimated 
from the integrated extragalactic contribution to 
dark matter annihilations, and compared with the 
measured atmospheric neutrino flux. Unfortu- 
nately, the extragalactic flux strongly depends on 
the shape of dark matter halos and their degree of 
dumpiness. A robust estimate is thus difficult to 
achieve. Although in Ref. [l| a relatively modest 
value of 2 X 10^ for the enhancement due to the 
dumpiness of DM was used, even values lower by 
a factor of 0(10) are possible. To be more con- 
servative, we use the diffuse photon flux from the 
smooth DM distribution in the halo of our Galaxy 
since: (i) Its normalization and distribution is bet- 
ter known (within a factor ~ 2); (ii) It is truly a 
lower limit for the DM annihilation flux Q . Sub- 
structure in our halo is expected to augment it up 
to orders of magnitude (see e.g. the parametric 
study 8] for our Galaxy or the study 9] for dwarf 
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galaxy satellites). Note that the contribution from 
the diffuse extragalactic photon background from 
DM annihilations further enhances the total DM 
emission. By neglecting both the substructure in 
our halo and the extragalactic contribution, we are 
being conservative. 

The differential flux of photons from dark matter 
annihilations is [1^ 

hnAE,^}) = —— — — / ds , 

where r{s,'ip) — (rg + — 2 Tq s cos V')"'" ip is 
the angle between the direction in the sky and the 
galactic center (GC), r© « 8.0 kpc is the solar dis- 
tance from the GC, and s the distance from the 
Sun along the line-of-sight (l.o.s.). In terms of 
galactic latitude b and longitude I, one has 

cos -0 = cos 6 cos L (2) 

Particle physics enters via the DM mass mx, the 
annihilation cross section ((Tann'y)j and the photon 
differential energy spectrum dN^/dE per annihi- 
lation. Concerning the DM halo profile, we adopt 
for the smooth DM mass density psm a Navarro- 
Frenk- White profile [13] 




where we choose pq = 0.3GeV/cm^ as the dark 
matter density at the solar distance from the GC, 
and a = 45 kpc as the characteristic scale below 
which the profile scales as r~-^. The galactic halo 
DM flux has a significant angular dependence, with 
possibly large fluxes from the galactic center re- 
gion. However, the DM profile in the inner re- 
gions of the Galaxy is highly uncertain. To be 



conservative, we shall only use the NFW profile 
for r > 1 kpc, a region where numerical simula- 
tions of DM halos have reached convergence and 
the results are robust [111, Of course, other 

choices for the profile are possible, but all of them 
agree in the range of distances considered here, dif- 
fering primarily in the central region of the halos. 
Since here we are focusing on the galactic diffuse 
emission rather than that from the GC, the resid- 
ual uncertainties which are introduced through the 
choice of profile (a factor ~ 2) are negligible for our 
discussion. 



III. THE DIFFUSE GAMMA RAY 
BACKGROUNDS 

The overall diffuse gamma-ray radiation can be 
qualitatively divided into a galactic and an ex- 
tragalactic contribution. Since the latter is not 
simply the isotropic part of the flux, the sepa- 
ration of these two components can be done at 
present only assuming a specific model for the 
production of secondaries by cosmic rays in the 
galactic disk and halo. (However, a measurement 
of the cosmological Compton-Getting effect that 
should be achievable for GLAST would provide a 
model- independent way to separate the two contri- 
butions A significant fraction of the quasi- 
isotropic component, especially in the GeV range, 
may be due to high-latitude galactic emission com- 
ing from processes in the magnetized halo of the 
Milky Way. For our purposes here, a detailed anal- 
ysis is not required, and thus we employ a fit of the 
galactic diffuse flux proposed in 14] and calibrated 
on EGRET data around the GeV [3], 



I^^iiE) = No{l, b) X 10-*^ (g^) cm-2s-isr-iGeV-i , (4) 



where the arguments are in degrees, -180° < / < 180° and -90° < 6 < 90°, 
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0.5, m>30° 



M (1 - ) \/l + (V35)Vl+[''/(l-l+0.022|!|)]2 ■ ' M - 



Vl + (V35)2yi+(b7T 

The EGRET collaboration derived the intensity of the extragalactic gamma-ray flux as 

-2.10±0.03 

0.451GcV j 



I^^{E) ^ (7.32 ± 0.34) X 10"*^ ( 777777^777 ) cm-^s-isr^^GeV"^ , (6) 

vahd from E ^ 10 MeV to E ^ 100 GeV. 
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The reanalysis of the data performed in [^j, based of cosmic rays, deduced an extragalactic spectrum 
on a revised model for the galactic propagation 
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FIG. 2: Bounds on (aannw) versus mx from diffuse 7 
rays (blue arrows) , atmospheric neutrino data [J] (ma- 
genta line) together with the expectation for a thermal 
relic (for s-wave annihilation), the KKT model and the 
unitary limit. See the text for details. 



FIG. 1: EGRET data for the diffuse extragalactic 
gamma ray flux, according to Q, and the fit of the 
original analysis in 3]. 

significantly lowered with respect to Eq. ([6|) at in- 
termediate energies, while closer to the original re- 
sult of Eq. ((6]) at the lowest and highest energy 
points. In Fig. [l] we show the points according to 
this reevaluation, together with the fit of Eq. ^ . 
To derive our constraint, we shall ask that the pho- 
ton flux from DM annihilations, integrated in each 
of the energy bins of Fig. [T] and in the whole en- 
ergy range covered by EGRET, remains below the 
sum of the upper limit for the extragalactic flux 
plus the galactic emission estimated according to 
the fit of Eq. (jl]). Note that since we compare 
the signal with the sum of the two contributions, 
the precise extragalactic fraction of the diffuse ra- 
diation is basically irrelevant. To be conservative, 
we shall compare the DM photon flux to the back- 
ground profiles along the curve I — 0, since the 
galactic background is maximum at this longitude 
(see Eq. g])). 



IV. GAMMA-RAY EMISSION FROM DM 
ANNIHILATION INTO NEUTRINOS 



dominant 2 — > 3 and 2^4 processes are the 
bremsstrahlung of an electroweak gauge boson 
that subsequently decays: XX vvZ^ve^W^ 
and XX vvj] . The branching ratio 

R — a{XX iJvZ) / a{XX vv) depends gener- 
ally only for ~ n?-^ on the details of the under- 
lying 2 — > 2 process. One can distinguish three dif- 
ferent regimes of this process: i) the Fermi regime 
mx < mz with 0{R) = [a2/{4:Tr)]^{mx/mz)'^, 
a) the perturbative electroweak regime 
mz < mx < a2/(47r)ln^(mx/mz)2 ~ lO^GeV 
where R grows from C'(a2/(47r)) to 0(0. 1), and 
in) the non-perturbative regime where large 
logarithms over-compensate the small electroweak 
coupling a2 0]- Here, we consider regime ii) and 
can use therefore standard perturbation theory 
for the evaluation of R. Numerical values of R are 
given in Tab. 1. 

The dominant source of photons are tt^ produced 
in quark jets from W and Z decays. The resulting 
differential photon energy spectrum dN^/dE has 
been simulated using HERWIG [Ts'l . 

The obtained bound from the EGRET Hmit is 
shown in Fig. [2] with arrows together with the 
limit from Ref. [l| using atmospheric neutrino 
data. The upper extreme of the arrow indicates 
the bound obtained by comparing the emissions 



By assumption, the DM particles X couple 
on tree-level only to neutrinos. Hence the only 
possible 2^2 annihilation process is XX — > Di' 
with an unspecified intermediate state that has 
negligible couplings to SM particles. Then the 



TABLE I: The branching ratio R = a{XX 
T/vZ) / a(XX uv) as function of mx- 



mx/GeV 


100 


300 


1000 


3000 


10'* 


R/% 


0.01 


0.02 


0.87 


1.9 


3.4 
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at the highest galactic latitudes (6 = 7r/2, I = 0), 
while the lower extreme is the bound coming from 
the inner Galaxy emission {b — 1/8, ^ = 0). 
The length of the arrow thus quantifies the im- 
provement due to our simple, angular-dependent 
analysis. Indicated are also the required value 
for a standard thermal relic with an annihilation 
cross section dominated by the s-wave contribu- 
tion, (CTannf) ~ 2.5 X 10~^^cm'^/s, the unitary limit 
(cannw) < 47r/(wTO^) for V — 300km/s, appropri- 
ate for the Milky way, and the constraints on the 
cosmological relativistic energy density from (iGj . 

V. DISCUSSION AND CONCLUSION 

In this paper we have shown that, even if dark- 
matter particles annihilate at tree-level only into 
neutrinos, diffuse gamma-ray data provide inter- 
esting constraints on their annihilation cross sec- 
tion because of electroweak bremsstrahlung. These 
bounds are comparable to the atmospheric neu- 
trino bound from Ref. [l[ in the mass range be- 
tween ^ 100 GeV and the onset of the stronger 
unitary bound around 10 TeV. Any appreciable 
branching ratio at tree level in electromagnetically 
interacting particles would lead to much stronger 
constraints from gamma-rays, but they are not 
as conservative as the bounds derived here or in 
Ref. [l|. A major improvement in the gamma-ray 



bound is expected from the GLAST satellite [13], 
to be launched by the beginning of 2008. In partic- 
ular, GLAST should resolve most of the diffuse flux 
of astrophysical origin, and map both the galac- 
tic and extragalactic diffuse emission with much 
higher accuracy, thereby improving the bound de- 
rived here. On the other hand, our results also 
suggest that the neutrino bound may be tightened 
as well by considering the DM annihilation in the 
galactic halo and taking into account the strong 
angular dependence on the halo signal [2l| . 

As a further application of our results, we note 
that the electroweak higher-order corrections dis- 
cussed here also contribute to increase the robust- 
ness of the bounds on strongly interacting dark 
matter from the Earth's heat flow in Ref. [19]. 
Above the TeV scale, electroweak bremsstrahlung 
put a lower bound of 0(1%) on the energy released 
in othcr-than-ncutrino channels, thus guaranteeing 
that an appreciable energy is release by annihila- 
tions in the interior of the Earth even for models 
with tree-level annihilations in neutrinos only. 
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